Linage and Outage Intervals in 

Transmission Systems 

Composed of Links 

By S. O. RICE 
(Manuscript received April 3, 19G3) 

This note is of the nature of an addendum to a recent -paper on satellite 
communication systems. It is concerned with the distribution and average 
durations of innages and outages occurring in transmission systems com- 
posed of a number of links. The links of such a composite system may be 
either in series, as in a radio relay system, or in parallel, as in a many- 
satellite stjstem. Several results regarding composite transmission systems, 
including some due to D. S. Palmer, are reviewed, restated, and extended. 

I. INTRODUCTION 

This note is in the nature of an addendum to a recent paper of mine 
on satellite communication systems. 1 It is concerned with the same 
general problem, namely the reliability of transmission systems com- 
posed of links which fail independently. Various published results are 
reviewed and extended. A large part of these results is due to D. S. 
Palmer, 2 whose excellent work was overlooked in my satellite paper. The 
approach given here differs somewhat from that used by Palmer. 

Incidentally, questions similar to those discussed here have also 
appeared in connection with coincidences in counting devices. 

The notation to be followed is illustrated in Fig. 1. Suppose that a link 
in a transmission system is always either in one or the other of two pos- 
sible states, state (a) or state (b). For example, if the link is a satellite, 
(a) may be taken as the state of being out of sight and (b) the state of 
being visible. Again, if the link is one of a series of links in tandem making 
up a transmission line, we may choose (a) to be the state of working 
order and (b) the state of breakdown. In a satellite system the links are 
in parallel and in the transmission line they are in series. Fig. 1 applies 
to both cases. 

The light portions of the top line in Fig. 1 represent the intervals dur- 
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Fig. 1 — Combination of k independent alternating sequences to form the 
resultant alternating sequence. 



ing which Link No. 1 is in state (a), and the heavy portions the intervals 
of state (b). Similarly, the second and third lines represent the state 
intervals of Links No. 2 and No. 3. This is a fc-link system with k = 3. 
The last line represents the system state intervals. In system state (ak) 
all k links are in state (a). In state (bk) at least one link is in state (b). 
State (ak) corresponds to the "intersection" of type (a) intervals and 
state (bk) to the "union" of type (b) intervals. 

For the satellite system, states (ak) and (bk) correspond to "out- 
age" and "innage," respectively. For the transmission line they cor- 
respond to "working order" and "breakdown." This reversal of interpre- 
tation for links in parallel and for links in series has been mentioned by 
Palmer. 2 

The problem is to find the distributions of the durations t ak and hit of 
states (ak) and (bk). The lengths t a , tj, of the intervals shown in Fig. 1 
are supposed to be independent random variables with given probability 
densities p a (t), Pb(t). Usually p a (t), Pb(t) will be the same for all k 
links, but in the more general case the densities associated with the ith 
link will be denoted by p* (t), Pb (t). The links are assumed to operate 
independently of each other. It is also assumed that the system has been 
operating long enough to reach statistical equilibrium. 

Although t a and h are independent, t a k and tbk need not be. An example 
is given just below equation (25) in Section IV. 

The results given here do not apply to the case where the pattern of 
intervals in two or more links shows periodicities. For example, if all 
type (a) intervals of Links No. 1 and No. 2 are of length 1 and all type 
(b) intervals are of length 3, then (depending on the relative phase) 
there may be no type (ak) intervals and just one infinitely long type 
(bk) interval. 

The distribution of t„k depends only on p a (t). It is obtained in Section 
II for general p a (t). The expected value tbk of tbk depends only on the 
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expected values t a , h and is given in Section III. At present there seems 
to be no practicable method, other than simulation on a high-speed 
computer, of obtaining the distribution of hk for general p a (t), Pb(t). 
For exponential p a (t) a method due to Palmer 2 and Takacs (outlined 
in Ref. 1) may be used, but even this is difficult unless p b (t) is also 
exponential. This method is developed in Section IV and illustrated in 
Section V. Sections VI and VII are concerned with the special case k = 2 
but general p a (t), Pb(t). Now the determination of the distribution of 
t b2 depends upon the solution of an integral equation. A vexing problem 
which I have been unable to solve is to show that when p a (t) is exponen- 
tial the integral equation leads to the same distribution as does setting 
k = 2 in the method of Section IV. 

I am indebted to John Riordan, David Slepian, and Lajos Takacs for 
helpful comments. 

II. THE DISTRIBUTION OF l a k 

It is convenient to set 

Fait) = l Pair) dr (1) 

J t 

AM - f F a (r)dr/t a (2) 

L - f rp a (r)dr = f F a {r) dt. (3) 

Here F a (t) is the probability that t a > t, t a is the expected value of /„ , 
and A a (t) is closely related to C. Palm's 3 "next-arrival" distribution. 
If a(s) is the Laplace transform of p a (t), i.e. 

a(s) = P e-'paiOdt (4) 

then 

J e Fa\t)dt- 

For the special case p a (t) = ae " 

Fait) = AM) = e~ u ', l = 1/a 

(6) 

a(s) = a/(a 4- s). 
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To interpret A a (t), consider Fig. 2, which shows a line in Fig. 1 cor- 
responding to a typical link. Choose a point t = x at random [this means 
that when the choice is from the very long interval (0,T) the chance 
that X falls between t, t -\- dt is dt/T], Let I be the distance to the end 
of the interval (which may be of either type) in which x falls. Then 2 ' 4,5 
A a (r) is the probability that / > t, given that x fell in an (a) interval. 

It should be noted that expression (2) for A a (t) holds even when suc- 
cessive t a 's and h's are correlated. This point is important in the proof 
of (7), since the intervals t ak , hk may be correlated. The only require- 
ment is that as T — * °° the distribution of the lengths of the (a) intervals 
in (0,2') approaches a definite distribution F a (t) possessing an average 
l a which is neither zero nor infinite. For emphasis we sketch a proof of 
(2) which is tailored to Fig. 2, in which, for the moment, /„ and k may 
be correlated. The chance that r < I < r + dr is the limit as T — > °o 
of the ratio 

[number of (a) intervals longer than t in (0,T)](dr) 
total length of (a) intervals in (0,T) 

In the limit this ratio approaches NF a (r)dT/Nl a , where N is the number 
of (a) intervals in (0,2'). Cancelling the N's and integrating t from t to 
oo then gives (2). 

To find the probability F„k(t) that /„* > t, suppose that all links are in 
state (a) at the randomly chosen time x. Since the links are independent, 
the chance that none has changed to state (6) by time x + I is [A a (2)] • 
Hence the function A ak (t) corresponding to the complete system is 
[A a (t)] h . This A a k(t) is related to F ak (t) by an equation obtained from 
(2) by replacing the subscripts ,l a" by "ah." Differentiation gives 

FAD = -l*£liMt 

dt ( 7 ) 

- (t uk /t a )k{A (t)\ k - l I<\(l) 

where l a k denotes the expected length of intervals of type (ah). 
Setting t = Oin (7) and using F ak (Q) = A a (0) = F a (0) = 1 leads to 

tak = ta/Ic. (8) 






-i t^r\ 



t=o t=x 

Fig. 2 — A„(l) is the chance that I > t. 
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When the individual links have probability densities p« (,) (0> P<> ' (0* 
i = 1, 2, • • •, k, the chance that the length of a type (ak) interval 
exceeds t is 

F ak (t) = -ld t flAa lii (t) (9) 

at i=i 

which implies 

cur = E (u^r 1 do) 

just as (7) implies (8). The A a 's and i a 's are related by equations cor- 
responding to (1), (2) and (3). These results are due to Palmer, 2 who 
obtains (7) by a different argument. 

III. THE EXPECTED LENGTH OF INTERVALS OF VARIOUS TYPES, INCLUD- 
ING TYPE (bk) 

The expected value hk of the length of the type (bk) intervals is 
related to t uk by the equation 

&* = (-' -iU* (id 



(±-1)1 

\P«* / 



where p llk [not to be confused with the probability density p uk (t)] is the 
chance that the random point x will fall in an interval of type (ok). 
This follows almost immediately from 

tbk/l,k = Pbk/Pak , (12) 

P„k + PM " I- (I 3 ) 

A careful discussion of the probability p uk has been given by Weiss. 6 

A relation similar to (11) also holds for t b , t a , and the probability p a 
that the random point x will fall in a type (a) interval. Solving for p a 
gives 

Pa = Uta + t b y l (14) 

and when this is combined with p ak = p« fc , which follows from the in- 
dependence of the links, (11) becomes 



3\* 

lb 



fc.-(^-i)U-[(i+J) -1 

Palmer's generalization of (15) can be written as 



£. (15) 



2272 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1963 

4i-[(IIp. ( V-iIU 

where lak is given by (10). 

Einhorn has given the instances k = 2 of (10) and (16). He also 
gives a generalization in which all k links are alike and attention is fixed 
on the average length t aT ,k' of the periods during which r or more of the 
links are in state (a).* He takes both p a (t) and pb(t) to be exponential, 
but his expression for t ar ,k appears to hold for general distributions. 
Thus, let state j be the state of the system in which exactly j of the k 
links are in state (a), and let 

p«M = Q)p«y-' d7) 

be the fraction of time the system spends in state j. Here f . j is a bi- 
nomial coefficient, p b = 1 — p a , and p a is given by (14). Then Ein- 
horn 's results may be stated as 

k 

tar.k = tar ljPaj.k/Par,k (18) 



1= 
k 



-1/3 \*-r 



• E({)WU) w /rr)(U rt (« 



tar.k" = tar J2p,tj,k/Par.k (19) 

where i ar = i a /r and !,„/' is the average length of the intervals during 
which j < r. Setting r = k in (18) and (19) gives (8) and (1">), re- 
spectively. 

To establish (18) note that, in the very long interval (0,T), the 

k 

amount of time the system spends in states for which j ^ r is T ?.y a 4.k . 

The number of periods in (0,T) during which j ^ r is equal (to within 
one) to the number of periods during which j < r, and both are equal 
to the number of times the system jumps from state r to state r — 1. 
As shown in the next paragraph, the number of these jumps is Tp art k/i ar , 
and (18) follows by division. 

Divide the links into two groups, Group I consisting of the first r 



* A similar problem has been considered in unpublished work by my colleague 
H. Coo, in winch account is also taken of repairs at periodic intervals. 
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links and Group II of the la.st k — r links. The fraction of time all k — r 
links in II are in state {b) is Pb~ r . The number of times all links in I are 
in state (a) is p/T/t aT . This is also the number of times Group I jumps 
from state r to state r — I. Since the two groups operate independently 
and no periodicities exist, we assume that pb~ T gives the fraction of 
these jumps occurring while all links in II are in state (6). Thus \p a r T/ 
iar]pb k ~ r is the number of jumps the complete system makes from state 
r to state r — 1 when a specified set of k — r links (namely the last 

A: — /) remain in state (b). Since the set may be chosen in I, _^ J 

ways, the complete number of jumps is Tp art k/t ar , as stated. 
Incidentally, it may be shown that 

FajAO = -l^^AjiDAt^t) (20) 

(U*r ! = j(t,)~' + <* -./H*.) - 1 

give the distribution and average length of the state j intervals. 

IV. THE DISTRIBUTION OF tbk 

In the first part of this section several auxiliary distributions are dis- 
cussed. They correspond to arbitrary p a (l), pi{t) and are more general 
than needed here, where ultimately p a (t) is required to be exponential. 
However, they are used in Section VI. 

Consider the probability Q aa {t) that x + t falls in a type (a) interval, 
given that the random point x falls in a type (a) interval. We have 

QaM = AM) ~ f Piait - t) ~A a (r) dr (21) 

in which A a (t) is the chance that the original (a) interval lasts beyond 
x + t and [— dAa{T)/dT] dr the chance that it ends in x + t, x + t + dr. 
The end of the original (a) interval marks the beginning of a type (b) 
interval, and Pba(t') is the probability that a type (a) interval exists at 
time t', given that a type (b) interval began at time 0. 

Let the Laplace transforms of p a (t), Pb(t) be a(s), 0(s). Then the 
transforms of F a (t), A a {t) are given by (5). Weiss, fa and Brooks and 
Diamantides, 8 have shown that the transform of Pba(t) is 

•"'[] - a(s)]0(s)/[\ - a(s)0(s)}. 

This result is also developed in my paper 1 in ignorance of the earlier work 
of Weiss. Since the integral in (21) represents a convolution, its trans- 
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form is the product of the transforms of dA a (t)/dt and Pb a (t). When the 
transform of Q aa (l) is computed from (21), it is found to be 

<f> a (s) = I e~ 8t Q aa (t)dt 



22) 

_ 1 [1 - a(m - Pis)] 
s sH a [l - a(s)0(s)] ' 

This result is given by Palmer" and, independently, by Brooks and Di- 
amantides. It is also given, together with a number of related results, 
by Cox (Ref. 5, Ch. 7). 

The argument in the two preceding paragraphs is concerned with type 
(a) intervals. It applies equally well to intervals of type (ak) when the 
(ak) and (bk) intervals are independent. When the links are alike, in place 
of Q ail (t) we have [Q aa (t)] k for the chance that a type (ak) interval exists 
at time x -f- t, given that one exists at the randomly chosen time x. In 
place of the probability densities p a (t), Pb(t) and their transforms a(s), 
jS(s) we have p a k(t), pbk(t) and their transforms a*(s), @k(s). Equation 
(22) goes into an expression for the Laplace transform of [Q aa (t)] k 

i _ [i -«.(*)iu -Awi . r ,-[gu(or * (23) 

8 sHah[l ~ a k .(s)(3 k (s)] J 

Since ctk(s) may be computed from (7) and Q aa (t) from (22), j8*(a) 
is the only unknown in (23). In principle, if not in practice, (23) may be 
solved for j9fc(s) and then pbk(t) obtained by inversion. 

It should be remembered that (23) is based on the assumption that the 
(ak) and (bk) intervals are independent. They are independent when 
p a (t) is exponential, since then p a k(t) is also exponential, and a knowledge 
of the lengths of the (ak) intervals tells us nothing about the (bk) in- 
tervals and vice versa. 

Thus when p a (t) = ae~ nt , so that a*(s) = ka[ka + s]~ , (23) reduces 
to 



1 



- f~ e-'[Q m (t)f dt (24) 



s + ka — kafik(s) h 

where Q aa (t) has the transform [s + a — afi(s)]~\ More generally, when 
p a (t) (t) = a i exp (—ait) and pb (t) is arbitrary the equation for jS*(a) 
becomes 



st ak 4 i - p k (s) J fcS 



'""'" = [ e- 1 flQ aa (i) (t)dt (25) 

•'n ;=i 
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where l ak = [J^o,-] -1 and Q aa lt) (t) has the transform [s + a, — ai/3 (,) (s)]~ l . 
This is essentially the result obtained by Palmer 2 and Takacs. (Takacs' 
version is given in my paper. 1 ) 

The only case in which the (ak) and (bk) intervals are obviously in- 
dependent seems to be that for exponential p a (l). On the other hand, the 
following example shows that successive (ak) and (bk) intervals may be 
correlated even though the (a) and (6) intervals for the individual links 
are not. Let /„ and h be uniformly distributed between 1.00, 1.01 
and 2.00, 2.01 respectively. Let /.• = 2. Then, given an (ak) interval of 
length tai = 0.f>, we can infer that the length of the following (bk) in- 
terval lies between 2.5 and 2.52. 

Hence, so far as the discussion given in this section goes, (23) is no 
more general than (24). The question now arises as to the form taken by 
Pbk(t) when p a (t) and pb(t) are arbitrary. Some information on this is 
given in Section VI for the case k = 2. 

V. DISCUSSION AND EXAMPLES 

When k tends to infinity, with p a (t), Pb(t) fixed but arbitrary, In- 
tends to zero and t ba to infinity. When t becomes small, (2) shows that 
A a (t) tends to 1 — t/l a and its fcth power to exp (—tk/t a ). It follows 
that when t is small and A; is large, the chance that the length of an (ak) 
interval exceeds t is 

Fak(t) w exp (-tftak). (26) 

It may be conjectured that when k — > co the chance Fbt(t) that 
Uk > t also tends to an exponential 

JW0-* exp (-*/&»). (27) 

Here t is supposed to be many times larger than ~t a and h . Indeed, con- 
sider a (bk) interval to be in progress and k to be large. Over all of the 
interval, except for a negligibly small fraction near the beginning, the 
process will be uncorrelated with the initial conditions, and the chance 
that the interval will end in (/, t + dt) is independent of t. This leads to 
the exponential form (27). When t b /l a « 1, k may have to be extremely 
large before (27) begins to hold. This is because the argument pictures a 
great deal of overlap of (b) type intervals. 

Next we turn to the case where the A - links are different and 



Pa (i) (t) = a.,-<r a '', Ph u> = b*r hi \ i = 1,2, ■■■ , k 

t a = (Li , lb — 0,- . 



(28) 
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This is one of the few cases in which the work may be carried forward 
to even a moderate degree. From Section II 

F a {i) (t) = A a (i) (t) = e~ ait 

tak = [2 Oi] _1 

^ (29) 

Fak(t) = exp [-t/tak] = exp [-«£a<] 

Pak(t) = E oj exp [— tJ2 «•] • 

These expressions for F ak (t) and p a k(t) also hold when pb (l) (t) is arbi- 
trary. From Section III 

kk = I 1 J tak , 

7* ' (30) 

Pak = II P« (<> , Pa U) = bi(Oi + 6,)" 1 . 

t = l 

The first step in using (25) is to compute Q aa ' (t) by inverting its 
Laplace transform. The result is given in the last line of Table I in Sec- 
tion VI and leads to 

fl a.'"(0 = II I" *' + of""""' ] = E c<e-*" (31) 



where in the general case the sum on j contains 2 terms. The right-hand 
member of (25) becomes J^ Cj(s + dy)~~\ and it follows that the La- 
place transform of Fu,(t) is equal to 

1 - k (s) _ 'tak . 

^E^ + ^r 1 (32) 

i 

The transform of Fbk(t) is thus a rational function of s. Its poles are 
at the zeros of ^2cj(s + dj)~\ and these zeros lie on the negative real 
axis between the points s = — dj, the rightmost of which is s = 0. The 
nth moment, tbk n , of hk is w!( — l) n ~ times the coefficient of s"~ in the 
power series expansion of (32). Hence for n > 

5? ■ *• [{-£f IzJ+dsr' - j }L ■ (33) 

When the links are alike, some of the d/s are equal, and 

[a , — (n+b)t-lk k /, \ it— n n 

l> + ae y* / k \ b a -„( B +mi 

b + a J "iWtHa) 4 
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Iii this case, the results are those used in Ref. 1, namely 
F*(t) = e~ tka , tak = 1/ka 
FhM = (1+p) w g exp [(a + b)zJl (34) 

/,-p m=0 2 m J (2m) 

t bk = [(1 + p) fc - l]/*o» p = a/6 
where// (2) = df k (z)/dz, and 2 , z x , • • • , Zk-i are the zeros of 

AGO -t(*)rrr. (35) 

„_o yi/ 2 + n, 

These zeros lie between the poles at 0, — 1, —2, • • • , — k. The first few 
terms in the power series for F b k(t) are given by (14) of Ref. 1. In present 
notation 

_ m />/ [(fc - i)o + bl&e 2 
JWfl = l - y\ + v 

(36) 
- [(k - 1)V + 4(fc - l)ab + 6 2 ] !]•+•••• 

The nth moment, w > 0, is 

In particular 

2(1+,)'- I; AV (38) 

akb iZi \n/ n 

a result given by Palmer for a = b = 1. 



When k = 2 and the links are alike, (32) gives 
1 - 2 (s) s + a + 26 



(s) s' 2 + (36 + a)s + 26 2 



(39) 



Results of the sort given here have occurred in the reliability studies 
of R. S. Dick 9 and others. 

VI. THE DISTRIBUTION OF 4* FOR k = 2 

Here we consider the determination of pbk(t) for the case of k = 2 
links when p a (l) and pb(t) are arbitrary. The following probabilities. 
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which are related to those mentioned in the first part of Section IV, will 
be needed: 

P{t) = chance that an (a) interval exists at time t, given that an 

(a) interval starts at time 0; 
Q(l) = chance that an (a) interval exists at time x -f- t, given that 
the random x of Fig. 2 falls in an (a) interval. Q(t) is equal 
to the Qaa(t) of Section IV; 
R(t) dt = chance that an (a) interval ends in /, / + dt, given that an 

(a) interval starts at time 0; and 
S(t) dt = chance that an (a) interval ends in x + /, x + t + dt, 
given that the random x falls in an (a) interval. 

Some information regarding these probabilities is summarized in Ta- 
ble I. Here a(s) and /3(s), the Laplace transforms of p a (t) and pb(l), are 
written for brevity as a and /3. The entries may be obtained by the 
methods indicated in Section IV. Closely related results have been given 
in unpublished work by my colleague H. E. Rowe. 

An expression for the distribution of the length of a (62) interval which 
consists of, say, three (6) intervals can be obtained by examining Fig. 3. 
The probability that the (62) interval shown in Fig. 3 has a length be- 
tween / and t + dt can be obtained by integrating 

Pb(t — x 2 ) dx 2 -R(x 2 — x 3 ) dx 3 -pb(x 3 ) dt 

(40) 
S(t — xi) dxi-pb(xi — Xi) d.Xi-P(x 4 ) 

over the permissible values of the x's, namely ^ x* ^ x 3 ^ x 2 ^ X\ ^ t. 
Starting with the probability Pb(t)Q{t) dt that a (62) interval consists 
of one (6) interval and is of length t, t + dt, one can write a series for 
pu 2 (t) dt. The first term is Pb(t)Q{t) dt, and the later terms are multiple 
integrals of the type obtained by integrating (40). An examination of 
the series shows that 



p b2 (t) = pb(t)Q(t) 4- / dxtSit - xi)q{t,xi) (41) 

•'o 

where q(t,x) satisfies the relation 
q(t,xi) = / dx 2 p b (t — a- 2 ) 

• Pb(xi)P(x 2 ) + J dx 3 R(x 2 — x 3 )q(xi,x 3 ) \. 

The series for pb 2 (t) obtained by repeated substitution of (42) in (41) 
is sometimes useful for small values of t. 



(42) 
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LIKE RANDOM POINT 
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I* X 4 — 



l Ii » 



Fig. 3 — Sketch illustrating quantities in expression (40). 
E( 1 nation (42) is of the form 

q(t,x) = g(t,x) + / K(t - x,x - y)q{x,y) d.y, (43) 

where g and K are known functions, and is an integral equation to deter- 
mine q(t,x) in the region ^ x ^ t. In this region the values of q(t,x) 
along the line x = Xi are expressed in terms of its values on the line 
t = xi . 

It appears difficult to solve (42) for general p a (t) and pb(t). Two 
special cases are discussed in Section VII. 

It should be possible to verify that when p a (t) = a exp (—at), (41) 
and (42) lead to the same result as does (24) with k = 2. However, I 
have been unable to do this for general Pb(t). The problem may be stated 
as follows: show that the /3 2 (s) defined by setting k = 2 and Q aa (t) = P(t) 
in (24) is the Laplace transform of 

Vbi(t) = p„(t)P(t) 4- a f P(t - x^qit^i) dxi (44) 

where P(t) has the transform [s + a — a/?(s)] _1 and q(t,xi) satisfies 

q(t,xi) = / dx 2 p b (t — x 2 ) 

■\p b (x 1 )P(x 2 ) + « / o " dx 3 P(x* - x»)q(x lt x a )\. 

When the two links have different statistics, the distribution of the 
lengths of (&2) intervals winch begin with a (b) interval of Link No. 1, 
say, may be expressed as an integral similar to (41). However, there are 
now two integral equations similar to (42) which must be solved simul- 
taneously. 
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VII. SPECIAL CASES FOR Pfc 2 (0 

Here two special cases are given in which the integral equation (42) 
may be solved. The second case has been used to study the problem 
mentioned in connection with (44). 

(*) Exponential p b (t) and general p a (t). When p b (t) = b exp (-bt), 
(42) shows that q(t,Xi) is of the form/(.*d) exp (-60 where 

/(.rx) = b f Xi f/.r 2C - 6< "-' 2) |~6P(.r 2 ) + f o ' R(x 2 - x 3 )f(x,) dxA . (45) 

This goes into 

F(«) = — \-r [b<p P (s) + <p K (s)F(s)] (46) 

s -\- b 

where F(s), (p P (s), and <pn(s) are the Laplace transforms of /(f), ^"(0, 
and 22(0- Similarly, the transform of (41) is 

fe(«) = W* + b) + <Ps(s + b)F(s + b). (47) 

From (40), (47), and Table I 

,,,/ x b\l - a(s)\ ( , R -, 

R(s s . _b_ sfr[l - a(s 4- b)} ( , 

P ' M) ' s + b (s + 6) 2 L[s 4- 26 - 26«(s + 6)] ' 

Inversion of 182(a) and [1 - /3 2 (s)]s _l now gives p b2 (t) and F b2 (t). 

When the two links have different exponential p b (t)'s and general 
p a (l)'n, the two simultaneous integral equations mentioned at the end 
of Section VI may be solved by a somewhat similar procedure. 

(ii) p a (t) exponential and p b (t) the sum of two exponentials. For 

p a (t) = ae-', p b (l) = r,c '"'' + c,e hil , cibr 1 + c4h x = 1 (50) 
(42) shows that q(t,x) is of the form 

c^'hix) + c,e- h %{x). 
Instead offj(x),j = 1,2, it is more convenient to deal with 

g ( x ) = p( x ) + a \ P(x - y)f } (y) dy 

having the Laplace transform Gj(s). When Laplace transforms are in- 
troduced, (42) goes. into two simultaneous equations for (li(s + b } ), 
(!-,(s + bo). Solving these leads to 
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ft(«) = CiGi(s + &,) + C2 G 2 (s + lh) 



-1 



C2B1 j- ciB 2 + 2acic 2 (s + bi + 62) 
£ii3 2 - a 2 dc 2 (s + 61 + 6 2 )- 2 ' 

R = 1 acj / v _ 1 (51) 

7 «P/-(s + bj) 8 + 26,- ' m ' s + a - a0(«) ' 

rf.) = * + C2 



s + 61 s + 6 2 
The equation for 0->(s) obtained from (24) in this case is 

[s + 2a - 2a i 9 2 (s)] -1 = / e~ st P\t) dt 
•'0 

= S— • / ^p( s — s'W(s') <**' 

2?rt ■'t 

(52) 
_ <pp(s)bib 2 . <pp(s — si)(si + 61) («i + bj) 

SiS, ' Si(Si — s 2 ) 

, y f (s — s 2 )(s 2 + &i)(«2 + ba) 

S 2 (S2 — Si) 

where the patli of integration L runs from — i<x> to +«*<» so that the 
singularities of <pr(s — s') lie on its right and those of <p(s') on its left. 
The integral has been evaluated by writing <p P (s) as 

, s (s + 6i)(« +h) 
<p,.{s) = 



s(s — Si)(s — s 2 ) ' 

(s — Si) (s — s 2 ) = « a + (a + 61 4- 6 2 )s + 6i6 2 + a(6i + 6 2 — c x — c 2 ) 

closing L by an infinite semicircle on the left and evaluating the residues 
at the poles s' = 0, 81 , s 2 . 

The task of verifying that (51) and (52) give the same value of fc(s) 
appears to be a lengthy one and has not been carried out. Several numeri- 
cal checks have been made and show no discrepancy. 
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